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ory for rigid rods. A phenomenological square-well potential is used to model these 
long-range interactions, and the model is used to compute phase separation and length 
fractionation as a function of well depth and rod concentration. The model closely 
captures experimental data points for isotropic/liquid crystalline phase coexistence of 
single-walled carbon nanotubes (SWNTs) in superacids. The model also predicts that 
the isotropic-biphasic boundary approaches zero as acid strength diminishes, with the 
possibility of coexistence of isotropic and liquid crystalline phases at very low concen- 
trations; this counterintuitive prediction is confirmed experimentally. Experimental 
deviations from classical theories for rodlike liquid crystals are explained in terms of 
polydispersity and the balance between short range repulsion and long range attrac- 
tions. The predictions of the model also hold practical import for applications of 
SWNT/superacid solutions, particularly in the processing of fibers and films from liq- 
uid crystalline SWNT/superacid mixtures. 

1 Introduction 

When suspended in a liquid, rodlike molecules and particles transition from an isotropic, dis- 
ordered state to an aligned, liquid-crystalline nematic state at sufficiently high concentration. 
This classical, peculiar behavior has generated numerous theories for predicting the isotropic- 
nematic transition and the equilibrium properties of the coexisting phases; typically, these 
theories treat the rodlike molecule or particle as an idealized Brownian rigid rod of high 
aspect ratio. These theories capture with some success the behavior of classical lyotropic 
liquid-crystalline materials such as liquid crystalline polymer solutions and rodlike biological 
particles in water [1]. However, previous theoretical studies fail to describe dispersions of 
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novel anisotropic nanomaterials such as single- walled carbon nanotubes (SWNTs). Here we 
extend the theoretical description of polydisperse Brownian rods to include the competition 
of short-range repulsion and long-range attraction that is known to yield rich behavior in 
colloidal phases. We focus on the coexistence between the isotropic and nematic phases. 

1.1 Isotropic-Nematic Phase Separation 

In 1949, Onsager developed a statistical mechanical theory for the rod orientation distribu- 
tion function in order to describe the phase behavior of idealized monodisperse rods inter- 
acting via a mean field excluded- volume potential in an athermal solvent p]. This theory is 
based on the competition between the orientational entropic contribution to the free energy 
and the excluded-volume potential. This competition gives rise to the isotropic-nematic 
phase transition that is characteristic of many rodlike molecules. As rods crowd beyond 
a critical concentration, alignment occurs spontaneously because orientational ordering is 
accompanied by a loss of interaction free energy which offsets the partial gain of free energy 
due to loss in orientational entropy. The predicted coexistence of an isotropic phase and 
an aligned phase is observed in experimental systems. Onsager's original formulation was 
numerically refined by Kayser and Raveche [3] and Lekkerkerker et al [3]. Later studies 
extended Onsager's approach to nonhomogeneous systems in order to resolve the profile of 
the interface between bulk isotropic and nematic regions El E] and finite-sized regions 
Pj. Onsager's approach was also modified to include attractive interactions via a simple 
order-parameter-based attractive potential by Khokhlov and Semenov [HI [ID] ; Khokhlov and 
Semenov used the Parsons scaling approximation to allow Onsager's second virial coefficient 
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truncation to remain valid up to the high concentrations that result from the attractive 
potential. 

The major competing theoretical framework for rigid rod phase behavior is Flory's lattice- 
based theory [TT], which accounts for excluded- volume interactions via packing effects and 
allows for attractive inter-particle interactions governed by the parameter x- Flory's predic- 
tions for the boundaries of the well-known biphasic chimney for an athermal solvent (x = 0) 
are slightly higher {7.89d/L, 11.57d/L) than the Onsager predictions {3.29d/L, A.19d/L), 
where d and L are the rod diameter and length, respectively. As the attractive force in- 
creases, the biphasic region broadens and the concentration of the coexisting aligned phase 
grows dramatically such that the coexistence nematic phase is a solid. (The Onsager-based 
results of Khokhlov [9l[T0] show a similar broadening effect.) The lattice-based Flory theory 
allows for great versatility in describing a variety of systems, including mixtures and various 
types of molecules. Flory and others used the lattice framework to extend models to semi- 
flexible rods [12], mixtures of rods and flexible coils [13], rods connected by flexible linkers 
[2] , rods with flexible side chains |15J , and, most relevant to the present study, polydisperse 
rods [16j. However, the myriad numerical approximations inherent in Flory's treatment have 
limited the quantitative applications of the theory to anisotropic nanoparticles [H], partic- 
ularly those that interact differently at long and short ranges. Moreover, the "attractive" 
interactions in Flory's theory simply indicate the enthalpy of mixing. In a poor solvent, the 
rods tend to "attract" one another simply because of the poor compatibility between rod 
and solvent. This attraction only acts on molecules immediately bordering the rod in the 
lattice, i.e., at short ranges. 

There have been an number of other efforts to describe attractive interactions between 
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rods. A number of studies have focused on depletion forces caused by non-adsorbing poly- 
mers in solution with the rods. Surve et al. modeled the effects of adsorption in a Flory-like 
framework [18] , while Bolhuis et al. used particle-based simulations to simulate spherocylin- 
ders with depletion forces caused by the addition of nonadsorbing polymer [12]. Borukhov et 
al. studied attractive depletion forces caused by cross-linkers in the context of stiff biopoly- 
mers in aligned phases [201 EI] ■ 

Particle-based simulations have also been used to simulate the phase behavior of rigid 
rods, usually described as ellipsoids, but these simulations are typically restricted to fi- 
nite aspect ratio [22] ■ The free energy of these systems are computed via Monte Carlo 
methods. Such simulations have been used to simulate the behavior of short aspect ratio 
spherocylinders interacting through attractive temperature-dependent square-well potentials 
[231 [2l] . Because of the shorter aspect ratio and the importance of temperature, these types 
of simulations are often applied to alkanes [22], far smaller than liquid- crystalline polymers 
or nanotubes; these simulations have also detected the existence of a nematic-liquid-vapor 
triple point for these shorter spherocylinders [2^ [27]. Particle-based simulations provide a 
useful means of evaluating approximations in Onsager-type frameworks; for instance, Sam- 
borski et al used Monte Carlo simulations as a point of comparison for various techniques 
of computing higher order virial coefficients for the Onsager theory |22j . 

1.2 Effects of Polydispersity 

In a system of rods, length polydispersity broadens the biphasic region, and long rods prefer- 
entially enter the nematic phase, as intuited by Onsager. Moreover, the coexisting isotropic 
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and nematic concentrations depend strongly on the initial (or "parent") concentration of 
rods. In a solution of monodisperse rods, these coexistence concentrations are indepen- 
dent of the parent concentration; in a solution of polydisperse rods, the increased degree 
of freedom in the length-dimension yields much richer phase behavior. Initial studies of 
polydispersity focused on long/short rod mixtures [Ij and thick/thin rod mixtures [7]. Sim- 
ple bidispersity in length or diameter can have important effects, including the formation of 
non-monotonic density profiles at the interface [7j. 

Three papers in 2003 rigorously analyzed the effects of polydispersity for continuous 
distributions of rigid rods. Speranza and Sollich computed the coexistence concentrations 
of a system of Onsager rods with continuous polydispersity in length [28l [29] . and they 
paid special attention to the effects of the longest rods in the distribution. Wensink and 
Vroege introduced an approximation for the excluded volume interaction and systematically 
computed the biphasic boundaries as a function of the initial concentration and the degree 
of polydispersity |30]. Notably, they found that, as polydispersity increases, the isotropic 
boundary moves to lower concentrations, and the nematic boundary shifts to dramatically 
higher ones, broadening the biphasic region. These three studies also identified specific 
scenarios where a triphasic (isotropic, short rod nematic, long rod nematic) equilibrium 
state can occur. Although the numerical analysis in [30] was limited to two forms of the 
parent population with varying polydispersity, the theory applies to systems described by 
any rod length distribution. 

These studies of polydisperse rod systems were limited to athermal solvents; here we aim 
to extend this work to include the effect of long-range attraction combined with short-range 
repulsion. 
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1.3 Experimental Motivation 

The present study is motivated by recent experimental results on the phase behavior of 
single-walled nanotubes (SWNTs) dispersed in superacids, such as fuming sulfuric acid and 
chlorosulfonic acid [211 E21 ESI [H]- Since lijima et a/'s observation of SWNTs in 1993 
[3^ , SWNTs have attracted great interest because of their unique electronic and mechanical 
properties [35]; much current research is devoted to the creation of neat SWNT macroscopic 
articles with properties comparable to the constituent SWNTs. However, the production 
of such macroscopic articles has been hindered by the difficulty of dispersing SWNTs as 
individual rods. 

The dispersion of SWNTs is difficult because they tend to bundle tightly due to strong 
van der Waals attraction to each other. Few solvents disperse pristine SWNTs as indi- 
vidual rods capable of forming a concentrated liquid-crystalline phase. SWNTs can form 
a liquid crystalline solution in water when wrapped or stabilized by acids containing hy- 
drophobic groups such as DNA or hyaluronic acid [2SIEZ| or when trapped in hydrogels [35] . 
Superacids, such as fuming sulfuric acid and chlorosulfonic acid, are the only fluids that 
spontaneously disperse SWNTs as individual rods in a concentrated liquid crystalline state 
[211 [221 HZ]- Superacids protonate the sidewalk of SWNTs and cause electrostatic repulsion 
and debundling [S^. These liquid crystalline dopes have been used for the production of 
macroscopic SWNT articles [2S1 [13110] • 

Surprisingly, in superacids weaker than (roughly) equimolar mixtures of sulfuric and 
chlorosulfonic acids, liquid-crystalline domains are in equilibrium with a dilute isotropic 
phase. These liquid- crystalline domains are long, threadhke structures, termed "SWNT 
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spaghetti," composed by a myriad of aligned SWNTs spaced by a few SWNT diameters. 
The SWNTs can translate along the axis of the strand, but have very limited mobility in 
the normal direction. This unusual phase behavior cannot be described by the Onsager 
or Flory theories for monodisperse rods because the SWNTs are polydisperse and because 
the theories cannot accurately describe the effects of varying acid strength on the inter-rod 
interactions. In fact, the broadening of the biphasic chimney occurs on the isotropic side 
rather than the nematic side, contrary to Flory's predictions. 

Modeling the phase behavior of such systems would involve both electrostatic forces due 
to protonation and inter-SWNT van der Waals attraction. The balance of electrostatic 
repulsive interactions and van der Waals attractive interactions for dilute polyelectrolyte 
colloids have been modeled using the Derjaguin, Landau, Verwey, and Overbeek (DLVO) 
theory [ITl 1^ . The DLVO approach was generalized for interactions between charged 
anisotropic macromolecules by Chapot et al |13] . The DLVO approach has been applied to 
a number of experiments on aqueous suspensions of charged rods such as V2O5 gels [Hj and 
b-FeOOH ^45j. DLVO and similar theories cannot be applied to rods in superacids because 
the assumption of dilute electrolytes inherent in the DLVO theory are inapplicable. 

We postulate that the key physical mechanism controlling the behavior of SWNTs in acids 
is the competition between long-range van der Waals attraction, which is weakly dependent 
on acid strength, and short-range electrostatic repulsion, which depends strongly on acid 
strength — as shown clearly by Raman spectroscopy of SWNTs in acids. We approximate 
this competition through a simple inter-rod repulsive-attractive potential, where strength 
of the attractive well varies with solvent quality. The rods interact through an attractive 
force that varies inversely with acid strength; this attractive force mimics van der Waals 
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interactions offset to varying degrees by electrostatic repulsion as a function of protonation. 



2 Formulation 

2.1 Free Energy 

Our formulation follows that of Wensink and Vroege and includes an additional potential to 
capture the effect of attractive interparticle interactions and solvent effects [30]. A solution 
of thin rigid rodlike molecules of diameter d and polydisperse lengths L is described by a 
rod distribution function A^(/). N{1) refers to the number of rods in the system with relative 
length /, where I is defined as L/Lq and Lq is an arbitrary reference length. For a volume V, 
N{l)/V is the local number density of rods of length I, and J N{l)dl gives the total number 
of rods in the system. 

The dimensionless Onsager excess free energy / for a volume V is written as 

/ = ~ / c(/)[lnc(/) - l]dl + / c{l)uj{l)dl 



VksT 

+ j j c{l)c{l')ll'[p(l,l') + \{l,l')]dldl\ (1) 

where the reference volume h is defined as irdLl/A, and the dimensionless rod distribution 
function c{l) is defined as bN{l)/V. The total rod number concentration is Cq, the zeroth 
moment of the distribution c(Z), where distribution moments are defined as 

Ck = J c{l)l^dl. (2) 

The overall concentration of the system can also be described by the rod volume fraction 
= {d/Lo)ci. The rod volume fraction is a more experimentally accessible quantity. 
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The entropic terms in the free energy description can be described as follows. The first 
term describes the ideal free energy of the polydisperse system. The second term contains 
the orientational entropy contribution, 



where ip{l,VL) is the normalized angular distribution function for rods with relative length 
I with orientation described by the unit vector u. This term is minimized by a uniform 
(isotropic) distribution of rod angles written as = I/Att. The third term contains the 
excluded volume interactions between rods. The average excluded volume contribution 
from rods of relative length I and I' is 



This term equals 1 for an isotropic state and is minimized by an aligned nematic state where 
the rods' alignment minimizes their excluded- volume interactions. The attractive/repulsive 
potential X{1, 1') models inter- rod attraction and repulsion, and is new. The phase behavior 
of idealized Onsager rods is governed by the balance between a;(/) and p{l, I') as a function 



We follow Wensink and Vroege in using the uniaxial Gaussian ansatz to approximate the 
angular distribution function ■?/^(/,u), written as 



such that ipiljU) is reduced to an unknown in only the /-dimension through the function a{l). 
The function a {I) is left as an unknown that describes the degree of alignment in a given 




(3) 





of c{l). 
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phase. Because of their higher excluded-volume cost, longer rods tend to be more ahgned. 
Wensink and Vroege write a;(/) ~ lnQ;(/) — 1 and then develop an asymptotic expansion of 
p{l, I') for the nematic phase in terms of a as 



p{l,l')r.^-{a{l)-^ + a{l')-'), (6) 

where terms involving higher powers of a{l)~^ are neglected. 

We designate variables as pertaining to the isotropic phase, nematic phase, and parent 
phase by using the superscripts (I), (N), and (0) for the theoretical description of a parent 
population c^'^^l) phase separating into c^^\l) and c^^\l). 

2.2 Equilibrium Conditions 

The degree of alignment for any nematic equilibrium state must minimize the free energy. 
In the present case, the degree of alignment is expressed through the unknown a{l), and the 
condition is written as 



Sa{l) 



dl' 



+ lc{l) I l'c{l')^-^dl\ (7) 



The excluded volume term in Wensink and Vroege's equation 7 is incorrect by a factor of 
two, while their equation 8 is correct. We argue below that the function A(Z, I') is a weak 
function of a{l) for the nematic phase such that the final term may be neglected. For an 
isotropic-nematic coexistence state to exist at equilibrium, two additional criteria must be 
satisfied. The two phases must have equal chemical potentials [p{l)/kBT — Sf/Sc{l)] as a 
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function of /, yielding the condition 

lnc(^)(Z) + 2lc^P = lnc(^)(/) + In [a{l)] - 1 + (8) 

and the excess chemical potential for the nematic phase jJiiP{l) is defined as 

/ii^HO = 2/ j l'c^''\l')\{l,l')dl' + 21 j c^''\l')l'p{l,l')dl\ (9) 

This criteria must be fulfilled for all /. The two phases must also have the same osmotic 
pressure [11 = 5// 5V] which yields the additional equation 

4'^ + (4'^)' = + j lc^''\l)l'c^''\l')W,l') + X{l,l')]dldl'. (10) 

Finally, conservation of mass must hold such that the two daughter phases in coexistence 
must add together to the parent population, i.e., 

c(o)(/) = (l-7)c(^)(/)+7c(^)(/), (11) 

where 7 is the volumetric fraction of the system that is nematic. Given a starting population 
c^^\l) and attractive potential A(/, /'), Eqs. ([71 [H [TOl [11]) comprise three functional equations 
and one scalar equation that may be solved for c^^\l), c^^\l), a{l), and 7. (Alternatively, 
the normalized parent population c^^^{l)/c^^ and 7 may be specified with the quantity Cq^"* 
as an unknown.) For convenience, we take the constant reference length Lq to be L^nean, 
the average length of the parent population. 

The Parsons scaling approach allows for use of the Onsager excluded volume at high 
concentrations fiO\. The difference in excluded volume is less than 5% for (j) < 0.2. We 
adapt the Parsons scaling approach for polydisperse systems to write the excluded volume 
term of the free energy as 

bFparson f ^^^^ ^ j^^^^ _ _ (13) 
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After some mathematical manipulation, the contributions to the chemical potential and 
osmotic pressure are then written as 

f c{l)c{l')ll'pil,l') 

>->-Parson ~ J ^-^ _ ■ (l-^j 

(t) ^ / ^" ~ "^^^^^p^^^ + {i-m) I ^^^^ 

The equation for equilibrium alignment becomes 



1/2 



Try V / 2a2(/) 



Lo\ (-ln(l-0(/)) 



l'c{l') 



._n-i/2 , 5A(/,/') 



(a(0-' + a(O"') + 



dr, (15) 



2.3 Attractive and Repulsive Interactions 



The novel contribution to the free energy expression is expressed in the quantity X{l,l'). We 
propose a particular phenomenological model inspired by the unusual "SWNT spaghetti" 
hquid crystaUine phase illustrated in Figure 1; the SWNTs in these strands are aligned along 
the strand with small distances between SWNT sidewalls. The SWNTs are no longer in 
due to the short-range electrostatic repulsion, but the attractive forces keep them within a 
few diameters of each another, such that each SWNT interacts with an effective "cage" of 
other SWNTs within the liquid crystalline phase. 

Our proposed model is as follows: A rod in an aligned state is subject to an attractive 
interaction with surrounding rods u{D), which represents the energy per unit length of the 
rod as a function of the distance D between the centers of the ahgned rods. This quantity 
scales linearly with I and is made dimensionless as U{D) = u{D)L^ean/ ^bT . The square- 
well potential U{D) is depicted in Figure 1 and is characterized by the four parameters C/max, 
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t^nin, Di, and D2. The discontinuities in the derivatives of the square- well potential are 
made continuous by using an arctan function to approximate the transitions of across 
Di and D2. 

As a first approximation, this attractive interaction is applied only to the dense nematic 
phase. This is motivated by experimental observations of the unusual liquid-crystalline 
phases in SWNT/superacid systems. The liquid crystalline phase is marked by ahgned 
SWNTs in long threadlike domains with acid molecules in-between, and the inter-SWNT 
distances are on the order of 10 — 30nm. This means that rods in the nematic phase are 
separated by average distances that are much smaller than the average distances between rods 
in the isotropic phase. This means that the rods in the nematic phase will spend substantially 
more time in the attractive well than rods in the isotropic phase; we approximate these 
experimental observations simply by applying the attractive well to the nematic phase but 
not the isotropic phase. The attractive interaction varies weakly with alignment in nematic 
phases. On average, the nematic phase will behave much like a columnar phase where the 
average inter- rod distance D can be calculated from the rod volume fraction in the 

nematic phase, as 



As the "test" rod feels the attractive potential, the rod's nearest neighbors in the densely 
packed nematic phase will "shield" the rod from others that are farther out such that 
the attractive interaction is minimal beyond those nearest neighbors. (Also, those near- 
est neighbors will be the first to enter the repulsive region.) Consequently, the quantity 
lc{l) J l'c{l')\{l,l')dl' is rewritten as lc{l)h.{l) where A(Z) = 6U{D). (A similar argument is 




(16) 
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used by Wensink and Vroege in the case of average excluded volume in the nematic phase.) 
The square- well potential can be parameterized by the dimensionless quantities t/max, t^min, 
Di/d, and D2/d, which may vary with solvent quality In addition, the repulsive region 
{D < Di) increases the effective excluded volume by a factor of {Di/d+l)/2. The repulsive 
portion of the square-well potential is applied to both isotropic and nematic phases. 

2.4 Shadow and Cloud Phases 

The phase behavior of polydisperse rod populations is complicated by the fact that the 
parent population c^^\l) influences the properties of the daughter populations c^^\l) and 
c^^\l). This stands in contrast to the monodisperse case where the coexistence values of 
Cq^ and Cq^^ are independent of the parent concentration Cq*^ . These concepts are illustrated 
in Fig. 2, which depicts the dependence of Cq ^ and on A thorough explanation of 
cloud and shadow phases may be found in the appendix. 

2.5 Parent Distribution and Solution Method 

For the SWNT samples of interest, the parent rod population is well-described by a simple 
Schultz distribution 



with normalization factor and is truncated at some minimum and maximum The 
polydispersity of any population is described by the parameter a, deflned as 



,(0) 



Nl'exp[-{z + l)l], 



(17) 




(18) 



which is related to the parameter z in the Schultz distribution as a 



2 




-1 
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In a typical purified batch of HiPco SWNTs (HPR batcli 152.2), we find that the parent 
distribution of SWNTs is well-approximated through the Schultz distribution as follows: 
d = Inm, Lmean = 348nm, Lmin = VJnm (the SWNT purification process eliminates SWNTs 
below a particular cutoff), I/max = 8/im, and a = 0.6 [IH]. We use this as an example parent 
population. 

In general, we explore numerically the phase behavior of a SWNT-dispersion as a function 
of the attractive potential well depth f/min, concentration 0q, and of the polydispersity a. 
For all studies, we set ?7max ^ |f^iin| since variations in t/max do not substantially change 
the phase behavior as long as f/max > 1? and we set D2/d 1 because nematic phases where 
D > D2 do not differ from the athermal case. 

The /-space is discretized into Q points in a fashion similar to that described by Wensink 
and Vroege, and Eqs. ([71 [HI \lQi [H]) are discretized into 3Q + 1 nonlinear equations and 
solved simultaneously using Newton's method. Solutions are traced in CQ°'*-space (or 7- 
space), a-space and f/min-space using arclength continuation in the parameter of interest 
[TU] . The following results and discussion are specific to the Schultz distribution but should 
be representative of the behavior of other parent populations. (A typical value for Q is 70. 
Note that the points need not be evenly spaced.) 
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3 Results & Discussion 



3.1 Cloud curves 

For the example population of rods described above with a = 0.6, the isotropic and nematic 
cloud curves are computed as a function of t/min for the case where Di/d = 1 (Fig. 3). 
f^iin = corresponds to polydisperse athermal rods; the isotropic cloud point is (j)^^^ = 
0.0049, and the nematic cloud point is (j)^^^ = 0.177. 

As attractive interactions grow, the nematic phase becomes more thermodynamically 
favorable and more concentrated. Because Di/d = 1, the aligned rods can pack tightly 
together in the liquid crystalline phase, so the nematic cloud curve moves rapidly to the 
maximum packing fraction as Umm decreases. Setting Di/d > 1 essentially sets a minimum 
for inter-rod spacing and a maximum Such a case is shown as a dotted line in Fig. 3 
for Di/d = 2.08, where the increase in (p^^'' is arrested by short-range repulsion. 

The nematic cloud curve in Fig. 3 for Di/rf = 1 is similar to the classic Flory theory 
[11] and the Khokhlov extension of the Onsager theory ^ because the biphasic region pri- 
marily broadens by increasing the concentration of the nematic cloud curve 0*-^^ toward the 
maximum packing fraction (0.907) as the solvent quality decreases. In all these previous 
studies, the nematic concentration rapidly increases to a solid-like state because short-range 
repulsion is not included in the model. The Parsons scaling approach causes the cloud curve 
to level off, similar to the results of Khokhlov et al. 

Most interestingly, the isotropic cloud curve diverges to very small values at a critical 
value of [/j^jjj = —0.581. Because logc*^'^)(/) ~ [(p^^^U^i^, the quantity log^'-'^-' will become 
more and more negative as the attractive force in the nematic phase increases. Therefore 
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in systems where attractive interactions are stronger than this critical value, the dilute, 
isotropic phase must always coexist with some liquid-crystalline phase. Even at very small 
concentrations, the system is still in the biphasic regime; thus, some rods will form an aligned 
phase at equilibrium due to the strong attraction. A similar result was found by Speranza 
and SoUich; for distributions with small amounts of extremely long rods, the cloud curve 
decreased to vanishingly small concentrations [28] . This holds great import for experimental 
studies of rod dispersions, as discussed below. 

Figures 4 and 5 show the influence of polydispersity a on the isotropic and nematic cloud 
curves for Di/d = 1. 

As polydispersity decreases, the isotropic cloud curves move down to lower well depths 
and to the right to higher rod concentration; this is to be expected since the effects of length 
fractionation decrease as polydispersity decreases. This also means that f/jnin decreases as 
a function of a; this holds practical import for experimental scenarios where an isotropic 
solution is desired. If a particular solvent/rod combination corresponds to a value of f/min 
that is less than f/^j^, then even dilute solutions will be in the biphasic regime. If the 
polydispersity of the solution can be decreased and/or the average length increased, then 
the value of f/^j^ will also decrease; thus, polydispersity can be used as a control to ensure 

that ?7min > ^min- 

As polydispersity decreases, the nematic cloud curve moves down to lower well depths and 
to the left to lower rod concentration. For low a, the nematic cloud curve is a much weaker 
function of acid strength than the isotropic cloud curve. Thus, for broad distributions (e.g., 
a = 0.6), the nematic cloud curve varies more rapidly than the isotropic cloud curve while 
for narrow distributions (e.g., a = 0.4) the opposite is true. Note that in the monodisperse 



case, the nematic cloud curve actually decreases with the attractive well becomes stronger; 
this differs from the polydisperse case because the attractive well amplifies the effects of 
polydispersity and broadens the biphasic region. 



3.2 Comparison with experiment 

In previous papers, we used a series of experiments to establish a phase diagram for SWNT/acid 
mixtures [3T| [T7] which is displayed in Figure 6. Various SWNT/acid mixtures at a given 
SWNT concentration (denoted by diamonds in Figure 6) are centrifuged into isotropic and 
nematic phases, and the concentration 0^^^ of the isotropic phase is measureco and marked as 
black circles. The concentration (p^^^ grows with increasing fractional charge (better solvent 
quality). In contrast with the predictions of Flory-like theories, the dilute isotropic phase in 
poor solvents coexists with a concentrated liquid- crystalline phase (0^^^ ~ 0.1) rather than 
with a sohd phase (i.e., c/)'-^-* > 0.9). 

In order to compare the theoretical results with experimental results for SWNTs in 
superacids, a relationship between U^-ain and acid strength must be posited. Superacid 
strength is typically measured in terms of base-specific Hammett acidity [17] . In the case 
of S WNT/superacid systems, acid strength is related to the fractional charge per carbon 
caused by protonation of the nanotube sidewall [21] ■ (This quantity varies linearly with the 
wavenumber shift, dG, of the G-peak in the SWNT Raman spectra.) The van der Waals 
forces should be roughly independent of acid strength. However, as the degree of protonation 
increases, the repulsion effect increases. Thus, f/min, which represents the sum of the van 



^Concentration is measured using UV-vis-NIR absorbance.[3T] 
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der Waals forces and the repulsive forces, becomes shallower as acid strength increases]^ 
We approximate the dependence of the repulsive force on the acid strength by positing a 
simple linear relationship between dG and the attractive well depth Umin by comparing the 
theoretical and experimental results for the isotropic concentration after phase separationjf 

Experimental measurements for the strongest acid used (chlorosulfonic acid) indicate that 
0*^^^ is approximately equal to the theoretical predictions for athermal (Onsager) rods; thus, 
we match a Raman shift of dG = 25cm^^ to f/min = 0. The theoretical predictions for (p^^'' 
for [/min = —0.957 match the values of (p^^'^ from centrifugation experiments on SWNTs in 
102% H2SO4 (i.e., 100% H2SO4 with excess 2 wt% dissolved SO3) which has a Raman shift 
of dG = ncm^^. These two data points set the linear relationship between f/min and dG. 

Figure 6 combines computational and experimental results in a full phase diagram for 
SWNT/acid mixtures. The model (red open circles) is remarkably successful at predicting 
the experimental values (j)^^^ (black circles) at low concentrations IV7\. Therefore, the simple 
square-well potential captures the chief aspects of the physical interactions between SWNTs 
in acids. 

Additional experiments were performed using the methods of Rai et al to test the model 

[21] . The concentration 0^^^ of the isotropic phase in 120% H2SO4 was measured as a function 
^ A similar effect is observed in DLVO theory for a fixed attractive force but a repulsive force that increases 

with increasing solvent quality. 

■^A number of factors will influence the relationship between acid strength and J/min, and these factors 

may vary from one SWNT sample to another. C/min is independent of SWNT length and polydispersity, 

but experiments indicate that the solubility of a HiPco SWNT sample in a given acid varies with diameter 

distribution, the frequency of defects in the SWNT sidcwalls, and the chemical effects of the oxidation and 

purification process on the SWNT sidewalls. 
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of the starting concentration cf)^ '. A comparison between theoretical and experimental data 

points is displayed in Figure 7. Again, the agreement between the two shows that the 

fundamental physics of inter-nanotube interactions are captured by the theory. 

Experimental measurements indicate a nematic cloud point of </)^^-' ^0.1 that is relatively 

constant with respect to acid strength; however, the theoretical results indicate a slightly 

constant value of (j)^^'^ = 0.115. Unlike the case of low-concentration phase separation, the 

theoretical and experimental values for the nematic cloud curve do not match. However, 

further analysis of the theoretical results reveals why the values differ. In Fig. 8, the 

theory predicts that 7 ~ 0.99 for a system where 0^"-' = 0.1, i.e., 99% of the system (by 

volume) is nematic. The remaining 1% is isotropic and comprised of relatively short rods, 

as shown in Fig. 8. In an experiment, these short-rod- dominated isotropic regions would be 

found primarily in the defects between liquid-crystalline domains, and these small isotropic 

regions would be nearly undetectable by optical microscopy, rheology, or dynamic scanning 

calorimetry. Thus, the experimental points are best taken as measurements of the transition 

to a system that is 99% liquid crystalline by volume, and we may take the theoretical 

prediction of the nematic cloud point as correct. The theory predicts that the starting 

concentration must be pushed up to (p^^^ = 0.115 to completely eliminate these isotropic 

regions. However, simulation results also indicate that if all of the rods with an aspect ratio 

under 40 are removed, then the nematic cloud curve decreases below 0'-°^ ^ 0.1,^ 

"^We do not believe there is nematic-nematic separation in our experimental system. Wensink and Vroege 
indicate that for lognormal distributions, there is a fairly small window of a and 0^"^ where this is possible. 
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3.3 Applications 

The computational results show that a liquid-crystalline phase will always be present for 
acids with dG < 20.15cm~^ (f^min < f^min' where ^7^;^ = —0.581, equivalent to an acid 
mixture with a 9:16 volume ratio of CISO3H to H2SO4), regardless of the starting SWNT 
concentration. This means that in order to get a dilute, isotropic suspension of individual 
SWNTs in superacid, the solution must first be centrifuged and phase separated in order 
to remove the isotropic phase from the nematic phase. A number of CNT characteriza- 
tion techniques (such as rheological characterization of CNT length [18]) require a dilute, 
isotropic suspension, and these counterintuitive computational results show that a simple 
low-concentration dispersion of SWNTs in 120% H2SO4 will not yield a dilute, isotropic so- 
lution. The isotropic cloud curve can be moved down to lower well depths and to the right 
to higher rod concentrations by eliminating the longest SWNTs in a given sample; this would 
be useful for applications requiring a dilute, isotropic solution and could be accomplished 
via a number of experimental techniques. 

Also, one of the chief applications of SWNT / superacid dispersions is the formation of 
aligned articles such as fibers and films [391 [U] j our computational results have immediate 
applications for the processing of such articles. The fiber spinning process works as fol- 
lows: High-concentration, liquid- crystalline SWNT / superacid dispersions are mixed, and 
extruded. The viscous dispersion becomes more aligned due to the tension on the fluid 
during extensional flow. Finally, the dispersion is coagulated in a non-solvent bath where 
the acid is removed, and the SWNT dispersion solidifies to form an aligned fiber. The best 
fibers are produced by dispersions that are highly aligned and fully hquid-crystalline, i.e.. 
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above the nematic cloud curve. Any isotropic regions in the dispersions form misahgned, 
weak regions in the final fiber; these defects decrease the electrical properties and are the 
primary points of fiber mechanical failure. (Note that defects that are predominantly com- 
posed of small nanotubes have already been observed for liquid-crystalline MWNT systems 

m-) 

With this application in mind, the theoretical framework outlined above can be used 
to locate the true nematic cloud curve in order to ensure that the dispersion is fully liquid 
crystalline. In this case, the cloud curve concentration may be too high for practical mixing 
and fiber spinning, so the theoretical framework can also be used to test the results of 
manipulating the starting SWNT length distribution in order to decrease the nematic cloud 
concentration. In this case, the elimination of SWNTs under 40nm will decrease the nematic 
cloud concentration below 0*^°-* ^0.1 such that SWNT/superacid dopes of (j)^^^ ^0.1 will be 
fully liquid crystallinej^ 

Also, this theoretical framework has implications for high-concentration rheology of 
SWNT dispersions. The dependence of viscosity on concentration is a well-known ex- 
perimental test for liquid-crystallinity; the viscosity increases, goes through a maximum, 
decreases due to the formation of a liquid-crystalline phase, and finally begins to increase 
again at high concentrations [33]. This dependence on concentration is blurred and weak- 
ened as the biphasic chimney is broadened because the formation of aligned phase occurs over 
''This may be accomplished by a numerical of experimental techniques. These include washing the 
initial SWNT sample with a weaker acid (such as 98% H2SO4) to eliminate the shortest SWNTs. Also, 
fractionating a sample in the biphasic regime in a weaker superacid (such as 102% H2SO4) and removing 
the isotropic phase can be used to eliminate the shortest SWNTs. 
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a much wider concentration range. SWNT/superacid dispersions with weaker acids or poly- 
disperse length distributions will show a weak dependence of viscosity on concentration and 
may not show the maximum/minimum signature of liquid-crystallinity. Decreasing poly- 
dispersity or strengthening the acid solvent will sharpen and amplify the difference between 
the maximum and minimum in viscosity. 

The utility of the theoretical treatment outlined in this study is evident from the unex- 
pected nature of the results. Optical microscopy experiments and DSC could not detect the 
isotropic regions in the SWNT/superacid dispersions used for fiber spinning at 0^°-' ~ 0.1, but 
the theory predicts that these regions are indeed present; such regions have adverse effects 
on as-spun fibers. Also, it is entirely intuitive to reason that a dilute solution of SWNTs in 
120% H2SO4 would be dispersed as isotropic individuals, but the theory predicts that even 
these low concentrations are within the biphasic region because of attractive interactions. 

We expect that this analysis will pave the way for the modeling of other anisotropic 
materials such as multi-walled carbon nanotubes (MWNTs) or inorganic nanorods. The 
phase behavior of CdSe nanorods [50] has been assessed through purely qualitative means, 
and an effective theoretical framework for understanding these materials is needed. The 
phase transitions for solutions of polydisperse functionalized MWNTs were qualitatively 
characterized by Song and Windle [511 [52j. A theoretical analysis of these systems may 
allow a more quantitative understanding of these dispersed nanomaterials. 
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3.4 Conclusions 

We have developed an extension to the Onsager theory for polydisperse rigid rods developed 
by Wensink and Vroege in order to capture the balance of long-range attractive and repulsive 
forces observed for anisotropic nanomaterials in solution. Our work is particularly motivated 
by the recent quantitative experimental data for SWNTs dispersed as individual rigid rods in 
superacids. Our results indicate excellent agreement between the theory and experimental 
data for predicting phase separation at low concentrations and for predicting the biphasic 
chimney's broadening on the isotropic side. The theoretical results also hold important 
and surprising implications for a variety of SWNT/superacid experiments and applications, 
including the understanding of SWNT/superacid rheology and the processing of fibers and 
films from liquid crystalline dopes. 
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Appendix 

Two useful but oft-misunderstood concepts for delineating the phase behavior of poly- 
disperse rod populations are the cloud phase and shadow phase. For example, at small 
concentrations, a parent population of polydisperse rods in an athermal solvent will re- 
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main isotropic with no phase separation. As the parent concentration Cg ' is increased, a 
point is reached where an infinitesimal volume of nematic phase, termed a "shadow phase," 
forms while the vast majority of the parent population remains in the "cloud" isotropic 
phase. Thus, the shadow nematic phase coexists with the cloud isotropic phase, and the 
corresponding value of c^^ represents the lowest value of c^^ at which coexistence is possi- 
ble. The converse is true on the high-concentration side of the phase diagram; the shadow 
isotropic phase coexists with the cloud nematic phase, and the corresponding value of Cq'^ 
represents the highest value of Cq°^ at which coexistence is possible. On the isotropic cloud 
curve, the quantity 7 (the fraction of the system that is nematic) will approach zero, while 
on the nematic cloud curve, 7 will approach one. 

The isotropic and nematic cloud phases are the ultimate boundaries of the biphasic region 
in a system of polydisperse rods. Parent populations within the biphasic region will phase 
separate into coexisting isotropic and nematic phases, and the values of Cq^ and Cq^^ will 
vary with These concepts are illustrated in Fig. ??, which depicts the dependence of 
Cq^ and Cq^'' on c^^ . Note that the daughter phases appear to remain in the biphasic regime 
which suggests that they are not an equilibrium solution; this paradox is resolved by the fact 
that the daughters have a different length distribution than the parent phase, such that the 
daughter phases lie outside the cloud curves for their specific distribution. 

Figure captions 

Figure 1: (I) The combination of short range electrostatic repulsive forces (shown in 
blue) and van der Waals attraction will result in an attractive well (shown in red) potential 
between the two particles. However, the exact form of these functions is unknown for complex 
solvent and solutes such as superacids and SWNTs. (II) A simple, phenomenological square 
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well potential U{D) is used to capture this balance of repulsive forces and attractive forces 
between SWNTs in various acids. As the acid quality decreases, the well becomes deeper. 
(Ill) D is the distance between neighboring SWNTs in the liquid-crystalline phase. This 
simplified potential is inspired by the unusual liquid-crystalline order seen in "spaghetti" 
in SWNT/superacid systems [33]. (IV) A schematic of "spaghetti" geometry depicts the 
threadlike, aligned nature of these unusual liquid-crystalline domains. 

Figure 2: Computational results for an athermal solvent illustrating the dependence of 
(p^^^ and (j)^^^ on (p^^^ for polydisperse distributions for or a = 0.6, Lmean = 348(i. Squares 
indicate shadow points, and diamonds indicate cloud points. This stands in stark contrast 
to the monodisperse case where (j)^^^ and (p^^^ do not vary. 

Figure 3: The isotropic and nematic cloud points are depicted as a function of the well 
depth f/inin for a = 0.6, Lmean = 348(i. When Di/d = 1, the nematic cloud point rapidly 
increases to (p^^^ = 1 as f/min decreases. For the case Di/d = 2.08, (p^^^ remains constant 
at 0.210. The Onsager predictions for monodisperse rods are depicted as triangles (A). 

Figure 4: Isotropic cloud curves (f/min vs. (p) as a function of a for Di/d =1. As a 
decreases, the isotropic cloud curves move down to lower well depths and to the right to 
higher rod concentration as length fractionation effects diminish. The critical value of ?7min 
drops with decreasing a. 

Figure 5: Nematic cloud curve for Di/d = 1 for various values of a. The cloud curves 
increase to = 1 as attractive interactions increase, but this effect diminishes with decreasing 
a. 

Figure 6: SWNT/superacid phase behavior as a function of SWNT volume fraction 
and acid strength (measured by dG). Black symbols denote experimental results [I?]. Red 
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symbols refer to theoretical predictions, and red open circles (o) designate the theoretical 
predictions for compared with experimental measurements of denoted by black 
circles (•). Black/red diamonds (o) indicate the initial system concentration 0*^°^ for both 
experiments and simulations of isotropic (I) - liquid crystalline (LC) phase separation. The 
open red triangles (A) is the Onsager predictions for (p^^'' for a system of monodisperse hard- 
rods. A red line (-) represents the theoretical predictions for the isotropic cloud curve. 
Simulations are for a = 0.6, L^ean = 348(i, Di/d = 2.08. 

Figure 7: SWNT/120 H2SO4 mixtures of varying 0^°-' were prepared and phase sep- 
arated; the isotropic concentration 0*-°'' was measured by UV-vis-nIR absorbance. These 
experimental measurements showed a close match with theoretical predictions (cr = 0.6, 
Lmean = 348(i) for phasc separation as a function of 0'-°^ 

Figure 8: (Left) Volumetric fraction of system that is nematic (denoted as 7) as a function 
of 0^°-' for for [/inin = 0, (T = 0.6, Lmean = 348(i, Di/d = 2.08. The system reaches 7 = 0.95 at 
0*-°^ ^ 0.1 (noted as triangle), but the concentration must increase to 0'-°^ = 0.115 in order 
to reach 7 = 1 (the nematic cloud point, noted as diamond). Experimental data indicates 
a cloud point of 0'-'''' 0.1 but cannot detect the small isotropic regions filled with short 
rods that persist up to 0^°-' = 0.115. (Right) Average length of rods in the isotropic phase 
(denoted as {Li)) as a function of 0'-°^ Near the isotropic cloud point, (L/) approaches that 
of the parent population, but (L/) decreases dramatically as 0*^°-' approaches the nematic 
cloud point. Thus, near the nematic cloud point, the rods in the isotropic phase tend to be 
the shortest rods in the parent distribution, as confirmed experimentally elsewhere [19]. 
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